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Numerical simulations indicate that several different glass transitions, for polydisperse systems in 
both two and three dimensions, are characterized by diverging correlation lengths. These correlations 
are described by Ising-like critical exponents, and are associated with diverging, Vogel-Fulcher- 
Tamann, structural relaxation times. Related simulations of thermalized hard-disks indicate that 
the curves of pressure versus packing fraction for different polydispersities exhibit a sequence of 
transition points, starting with a liquid-hexatic transition for the monodisperse case, and crossing 
over with increasing polydispersity to glassy, Ising-like critical points. I propose to explain these 
observations by assuming that glass-forming materials contain twofold degenerate, topologically 
ordered clusters of particles, similar to the two-state systems that have been invoked to explain 
other glassy phenomena. This paper starts with a brief statistical derivation of the thermodynamics 
of thermalized, hard-core particles. It then discusses how a two-state, Ising-like model can be 
described within that framework in terms of a small number of statistically relevant, internal state 
variables. The resulting theory is remarkably consistent with the simulation data, but fails near the 
transition points where Ising symmetry must cross over to rotational symmetry. I also suggest a 
possible relation between Ising correlation lengths and the Vogel-Fulcher-Tamann formula. 



I. INTRODUCTION 

Among the deepest outstanding questions in materials 
theory is whether glass forming systems undergo truly 
thermodynamic phase transitions, and whether, if so, 
those transitions are characterized by diverging length 
scales. I argue here that recent numerical simulations, es- 
pecially those by Tanaka and coworkers 0, 0] , imply that 
this thermodynamic picture is correct; and I propose an 
Ising-like scenario within which glassy phase transitions 
and critical correlation lengths may be understood. 

Many current theories of amorphous materials assume 
that glassy slowing down is fundamentally a kinetic phe- 
nomenon that may, or may not, culminate at a sharply 
defined temperature or density. For example, facilitated 
kinetic Ising models have no nontrivial thermodynamic 
properties. [3[ Similarly, the random first-order transition 
theory (RFOT) 0, Q starts by computing a thermally 
activated relaxation rate that implies a diverging length 
scale. My own excitation-chain (XC) theory Q originally 
was introduced as a way of avoiding the weaknesses I per- 
ceived in RFOT Q; but the XC theory uses an analogous, 
intrinsically dynamic strategy, and produces a similarly 
diverging length scale (with apparently incorrect criti- 
cal exponents). If glassy systems do, in fact, undergo 
thermodynamic phase transitions, then that equilibrium 
behavior ought to be the basis for a theory of relaxation 
rates, and not the other way around. Theories of relax- 
ation mechanisms ought to start with a theory of the 
equilibrium state - if such a theory exists. 

The strongest evidence in favor of glassy, thermo- 
dynamic phase transitions is the recent discovery that 
several different glassy systems exhibit growing bond- 
orientational order, with correlation lengths that diverge 
at critical temperatures or packing fractions. In [l|, 
Tanaka and coworkers describe molecular dynamics sim- 
ulations of a variety of two and three dimensional sys- 



tems, including polydisperse hard disks and spheres. A 
related analysis by Mosayebi et al Q, for a three di- 
mensional, bidisperse, Lennard- Jones system, produces 
results that are consistent with those of [l| , and provides 
additional mathematical rigor and physical insight. This 
length scale may be the same as the one deduced indi- 
rectly in numerical simulations by Karmakar et al. @. 

The picture that emerges in [l[ and Q is that the di- 
verging, equilibrium correlation lengths £ are described 
by Ising-like critical exponents. Specifically, £ is pro- 
portional to t~ v with v = 2/d, where d is the spatial 
dimensionality, and t is a dimensionless measure of the 
distance from a critical point. For thermally controlled 
systems, t = (T — T c )/T c , where T c is the critical temper- 
ature. For hard-core particles, the temperature T is re- 
placed by the inverse of the packing fraction 0. For both 
d = 2 and 3, the structural relaxation time r a is found 
to be consistent with a Vogel-Fulcher-Tamann (VFT) re- 
lation, log(T Q ) ~ t , which may be further evidence of 
Ising-like universality. 

More evidence bearing on the existence of glassy phase 
transitions is summarized in Fig. 2 of 0] , where Kawasaki 
and Tanaka have plotted the pressure p as a function 
of cj> for simulated hard disks at a sequence of increas- 
ing percentage polydispersities A. As expected, the 
monodisperse system at A = % exhibits an appar- 
ently sharp transition between liquid and hexatic phases 
at <f> = 0.69. With increasing A, the transition points 
on the p(<f>) curves move to larger p's and $>'s, and be- 
come less and less distinct. They are invisible in the 
pressure data above A = 9%, which is the value of the 
polydispersity for which Tanaka et al. [l[ report a bond- 
orientational correlation length that extrapolates to in- 
finity at <f> = 0.787. The important point for present pur- 
poses is that the sequence of pressure curves in Fig. 2 of 
[2| appears to indicate a smooth crossover from a liquid- 
hexatic transition at A = % to Ising-like critical points 
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for A > 9 % - a qualitative change of universality class. 

In thermodynamic equilibrium, hard-core systems are 
characterized entirely by steric constraints and configu- 
rational entropy. They are the antitheses of the solv- 
able models of systems with infinitely long-ranged inter- 
actions that sometimes are used as starting points for 
glass theories, e.g. RFOT. In hard-core systems, there is 
no stored elastic energy, and therefore there can be no 
such thing as the vibrational entropy sometimes postu- 
lated to be responsible for crystalline ordering in such 
materials. Nor, for classical systems of this kind, can the 
entropy associated with kinetic degrees of freedom play 
any role in determining equilibrium structures; the ki- 
netic energy simply factors out of the partition function. 
Thus, hard-core simulations, and analogous experiments 
using hard-core colloidal particles, present a particularly 
clean theoretical challenge. In this paper, I consider only 
such systems. Accordingly, in Sec. HH I present a brief 
statistical derivation of the thermodynamics of thermal- 
ized, hard-core particles, which also serves to emphasize 
the role of properly chosen internal state variables. 

In the main part of this paper, I use this thermody- 
namic framework to develop an Ising-like model of disor- 
dered, hard-core materials. The basic ingredient of this 
model is a population of twofold degenerate, topologi- 
cally oriented, clusters of particles, similar to the two- 
state systems that have been invoked to explain other 
glassy phenomena. [13, [ll| In Sec. IIII1 I propose a ratio- 
nale for this model, and then, in Sees. IIVI and [V] show 
that it predicts both Ising-like, bond-orientational cor- 
relations, and a sequence of critical ordering transitions 
along the p(<p) curves. In Sec. I VII I suggest a possible 
relation between the Ising correlation lengths and the ob- 
served Vogel-Fulcher-Tamann behavior of the structural 
relaxation rates, thus tentatively adding support for the 
two-state hypothesis. 

In Sec. IIII1 however, I also argue that Ising symme- 
try must cross over to full rotational symmetry when the 
correlations become sufficiently long-ranged near critical 
points. As a result, the two-state theory must fail in the 
immediate vicinities of ordering transitions. Neverthe- 
less, this theory accounts remarkably well for the simu- 
lation results, and may therefore provide a unified point 
of view for understanding a wide range of glassy behav- 
iors. I conclude in Sec. IVIII with some remarks about 
that conjecture. 



II. STATISTICAL THERMODYNAMICS 

In preparation for developing a model of thermalized, 
hard-core, glass forming particles, we need to understand 
the statistical thermodynamics of such systems. The fol- 
lowing analysis is based on [l2[ and is almost, but not 
exactly, identical to that presented in [l]|. 

In the absence of interaction energies, the only exten- 
sive quantity available for describing a hard-core system 
is its volume V. (Throughout this paper, the term "vol- 



ume" means either three dimensional volume or two di- 
mensional area.) V must be a function of the entropy S 
plus a small set of internal state variables that govern the 
responses of the system to external forces. Denote these 
variables by A Q , a = 1 , ... n or, equivalently, by the 
set {A}. As discussed in [12|, the A a must be extensive 
quantities, or spatial averages of such quantities, each 
carrying its own entropy. In the thermodynamic limit 
of very large systems, the total entropy, say S(V, {A}), 
must approach its equilibrium value when the A Q ap- 
proach their own equilibrium values. I discuss the choice 
of these variables in Section HVl but, first, I consider only 
the general structure of the theory. 

For simplicity, assume that the kinetic degrees of free- 
dom of the particles plus the degrees of freedom of the 
heat bath in which the system is immersed constitute a 
single thermal reservoir at temperature 9 = ksT '. Denote 
the energy of this reservoir by Ur. Because the hard-core 
configurational degrees of freedom carry no potential en- 
ergy, Ur is the total energy of the system. Therefore, the 
first law of thermodynamics is simply Ur — — p V. 

The total entropy of this system, S(V, {A}), is the sum 
of the configurational entropy, say Sc, and the entropy 
of the reservoir, say Sr. Sc{V, {A}) is a constrained 
entropy computed by counting the number of configu- 
rations with fixed values of V and {A}. Conversely, 
V = V(Sc, {A}). In analogy to the notation of Ed- 
wards and coworkers [H|-[l6j], define the compactivity X 
by writing 



X 



dV 



(2.1) 
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Thus, the first law becomes 
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The second law, Sc + Sr > 0, is best written by using 
Eq. ([2.2[) to eliminate Sc, with the result 
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A a > 0. 



(2.3) 



As usual, the requirement that this inequality be satis- 
fied for arbitrary variations of external conditions implies 
that each of the terms on the right-hand side of Eq. (|2.3j) 
be separately non-negative. In particular, the first term 
is non-negative if 
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which means that, in the equilibrium limit, X — > 9/p. 

To interpret the remaining terms on the right-hand 
side of Eq. (|2.3[) . remember that, in a statistical sense, 
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the internal variables {A} describe only a fraction of the 
number of configurational degrees of freedom of the sys- 
tem. To account for the remaining internal degrees of 
freedom, write 



and 



Sc(V,{A}) -5({A}) + 5i 



V(Sc,{M) = V({A}) + V x (Sx) 



(2.5) 



= V({A}) + V 1 [Sc-S({A})\ > (2.6) 

where 5({A}) is the entropy associated with the internal 
state variables; V({A}) is the corresponding volume; and 
Si and V\ are, respectively, the entropy and volume as- 
sociated with all the degrees of freedom other than the 
{A}. Then, for each a, 



dV 



dV 



dA a dA a 
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dA a dA c 



F({A}), 



where 



J({A}) = V({A})-I5({A}). 



(2.7) 



(2.8) 



Inserting this result into Eq. (|2.3|) . satisfying the second- 
law inequality separately for each term in the sum over 
a, and taking the equilibrium limit, we find that 



A„ oc 



d 



^({A})^0. 



(2.9) 



Thus, F{{A]) is a "free volume," analogous to a free 
energy, whose minimum in the space of variables {A} 
locates the equilibrium state of the system. The values 
of the {A} are determined by the equations 



P dS/dA a 



9 dV/dA a 



(2.10) 



Note that Eg. (12.101) is what we would have found had 
we simply maximized the entropy S for a fixed volume 
V, and used X = 9/p as a Lagrange multiplier. The 
preceding derivation is more general in the important 
sense that it includes the residual quantities 5*1 and V±, 
which will play a role in the following analysis. 



III. TWO-STATE SYSTEMS 

The first step in constructing a model of hard-core par- 
ticles within the framework outlined in Sec. HI] must be 
to choose the internal state variables A a . In doing this, 
we must decide which of the degrees of freedom of the 
system as a whole are statistically relevant, and there- 
fore should be included among the A Q , and which can be 
included implicitly in the residual quantities Si and V\ 
defined in Eqs. (l2.5p and 



In some way, one or more of these variables must de- 
scribe the bond-orientational order observed by Tanaka 
et al..[l] Tanaka's principal innovation has been to look 
for spatial correlations, not between particle positions 
per se, but between the positions of particles in topologi- 
cally similar environments. Specifically, in what I believe 
to be the most definitive of their two-dimensional sim- 
ulations [l[, Tanaka et al. measured the time-averaged, 
complex, hexatic order parameter $e as a function of 
position and packing fraction 4>; and they computed the 
two-point correlation (^>g(r) ^1(0)} as a function of the 
separation r. From the latter quantity, they computed 
the correlation length £(</>), and found that it scaled as 
described in Sec. Has a function of t = (4> c — cf))/(j) Cl where 
the critical packing fraction <f> c has replaced T~ . Simi- 
lar results were obtained in three dimensions, where the 
relevant topological order parameter was found to be the 
degree of hexagonal-close-packed (as opposed to icosahe- 
dral) order. (See also fl7j|.) 

It might seem that a topological order parameter such 
as <&6 would necessarily be one of the A a . Indeed, 
is frequently used as an argument of a Landau free en- 
ergy, from which equilibrium states of two-dimensional 
systems are determined by a variational procedure for- 
mally identical to Eqs.(2.8) and (.2.10). (For example, 
see [Is].) However, we need to look more closely at such 
models before trying to use them to describe glass for- 
mation. 

In the monodisperse limit, Tanaka's hard disks un- 
doubtedly undergo the liquid-hexatic transition that has 
been studied intensely ever since the pioneering numer- 
ical simulations of Alder and Wainwright . [l9| The stan- 
dard description of such transitions in microscopically 
uniform materials is the two-dimensional melting the- 
ory of Kosterlitz, Thouless, Halperin, Nelson and Young 
[lO-iH; but this "KTHNY" theory may not be what 
we need to describe the most important properties of 
polydisperse, glass-forming liquids. KTHNY describes 
the melting of an hexatically ordered phase as a process 
in which a dilute population of disclination pairs under- 
goes a thermally induced unbinding transition, thereby 
destroying long-range oricntational order in a distinctly 
non-Ising manner. In contrast, a glass-forming liquid, 
well away from a KTHNY transition, cannot naturally 
be described by a population of disclinations. Even if it 
were possible to do so in some formal way, we know that 
the KTHNY analysis fails when that population becomes 
too dense, as must happen in the liquid phase. Thus, it 
should be more productive to construct a theory in which 
topological order emerges from within a liquidlike state, 
instead of, as in KTHNY, starting from a state with in- 
finitely long range order and asking how it melts. 

Accordingly, I propose that the fluctuating liquid state 
of a glass-forming material be visualized as one in which 
topologically ordered clusters of particles appear and dis- 
appear in a background of disordered, fluidlike particles. 
These ordered clusters may be favored by steric (or en- 
ergetic) interactions; in the liquid phase, they are disfa- 
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vored by the entropy of the system as a whole. As the 
pressure is increased, they come closer together, and the 
steric forces make it favorable for them to be aligned with 
each other. Thus, topological order grows with pressure. 

My main hypothesis is that these clusters are statis- 
tically most likely to be two-state systems. The glass 
literature contains many references to such systems. 
For example, in 1972, Anderson, Halperin and Varma 
[lOj based their theory of low-temperature anomalies in 
glasses on the hypothesis that "in any glass system there 
should be a certain number of atoms (or groups of atoms) 
which can sit more or less equally well in two equilib- 
rium positions." More recently, my colleagues and I have 
used a similar argument to justify the model of two-state, 
shear transformation zones (STZ's) that we have used in 
theories of amorphous plasticity. [ll| Twofold symmetry 
is especially important for present purposes because it 
is the Ising symmetry, and thus is consistent with the 
observed Ising-like critical exponents. Moreover, the fol- 
lowing argument in favor of twofold symmetry is inde- 
pendent of the specific nature of the ordering, or even 
the dimensionality, and thus may lead to the kind of uni- 
versality that seems to be emerging in the computational 
experiments. 

To see how two-state systems might naturally occur 
in theories of disordered materials, note that a spatially 
varying order parameter such as ^6 should be defined as 
an average over some coarse graining length scale. If that 
scale is too small, say, only one or two particle spacings, 
then ^6 w iH be large in some places and small in others; 
but changing the local orientation at an hexatic site, i.e. 
changing the phase of almost certainly increases the 
local volume (or energy), so that the system is rigid at 
particle-sized length scales. At the other extreme, if the 
coarse-graining scale is very large, and if we are not too 
close to an orientational ordering transition, then many 
different orientations will be degenerate in the sense of 
having the same volumes or energies, but the averaged 
magnitude of w iU be too small to provide useful in- 
formation. I propose, therefore, that the natural coarse- 
graining scale - the one that provides the statistically 
most relevant information - is the one for which the or- 
dered clusters are most likely to be twofold degenerate, 
simply because two is the smallest integer greater than 
one. 

This coarse-graining argument immediately points to a 
limitation of the theory. As an Ising system approaches a 
critical point, the correlations become long ranged, and a 
renormalization-group analysis like that used by KTHNY 
requires that we coarse-grain on increasingly large length 
scales. For hard disks in two dimensions, we eventu- 
ally restore circular ( "xy" ) symmetry [2l[ , and cross over 
into a regime where the KTHNY analysis again becomes 
valid. As a result, even for a polydisperse system, there 
must be a region near a critical ordering transition where 
the correlation length diverges according to the KTHNY 
prediction. This crossover region may be unobservably 
small as a function of packing fraction for large polydis- 



persities; but I think it must be there in principle. Con- 
versely, the crossover may also occur for a monodisperse 
system, because the liquid phase is intrinsically disor- 
dered away from criticality. 

The fundamental question regarding the two-state hy- 
pothesis is whether it can be derived systematically from 
a well defined description of a many-body system. I see 
no reason why such a derivation should not be possi- 
ble. Perhaps the many-body strategy presented recently 
by Yaida 25], which also concludes that glassy systems 
belong to the Ising universality class, is a step in this di- 
rection. However, in the next several Sections, I simply 
take the two-state model literally, and examine how its 
predictions compare with the simulation data. 



IV. ISING-LIKE MODEL 

A. Binary Clusters 

To describe the two-state picture mathematically, let 
N + and iV_ be extensive variables denoting the numbers 
of, say, "binary clusters" oriented in + and — directions 
with respect to some direction in space. Degeneracy re- 
quires that, when a cluster switches between + and — 
orientations, it continues to make the same contribution, 
say v*, to the volume V introduced in Eq. (|2.6l) . Just as 
in the STZ theory, the actual orientations denoted by ± 
need not be specified initially. In the STZ case, we usu- 
ally interpret the N± to be the numbers of zones whose 
orientations are more nearly parallel or antiparallel to an 
applied stress. In contrast to the STZ's (or the disclina- 
tions), there is no reason why the population of binary 
clusters should be dilute. A large ordered region at high 
compression may consist almost entirely of aligned clus- 
ters, whose specific orientation may be the result of an 
accidental anisotropy or a spontaneously broken symme- 
try. 

The way in which orientational order is propagated 
between neighboring positions in this model is via an 
Ising-like interaction, in which neighboring clusters make 
smaller contributions to the volume if their orientations 
are aligned than if they are opposite to each other. An- 
other way of thinking about this is that the neighboring 
clusters break each other's orientational degeneracies in 
such a way as to increase the probability of their align- 
ment. This steric effect, the analog of an Ising exchange 
coupling, is the reason why orientational order increases 
in response to increasing pressure. 

To make direct comparisons with the functions p((f>), 
we need one more internal variable to describe how the 
system as a whole expands and contracts in response to 
changing pressure, and how that behavior couples to the 
internal state of topological order. For this purpose, it 
is useful to introduce a population of, say, Nq "voids" 
occupying volumes vq . In order to play a role comparable 
the the N±, No must be a collective variable describing 
a property of groups of particles comparable in size to 
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the binary clusters, and vq must be a volume associated 
with more than just one, particle-size void. With the 
extra degree of freedom described by No, the model can 
make a transition with increasing pressure from dilute, 
liquidlike states with large populations of voids, to dense 
ordered states in which the voids disappear. Tanaka et al. 
show such voids in their Voronoi tiling patterns for two 
dimensional systems (see Fig. 2 in [![), and assert that 
these voids play a role in limiting the extent of hexatic 
correlations. That happens here as well. 



B. Volume 

In a first, mean- field statement of this model, the vol- 
ume V defined in Eq. (|2.6j) is 



y at N * v * + Nq Vq _ 



J 



2(N* + N ) 



(N 



Nt 



(4.1) 



where N* = N + + A_ is the total number of binary 
clusters, and is therefore proportional to the extensive 
number of statistically relevant, orientational degrees of 
freedom associated with the partial volume V and the 
partial entropy S. The pairwise interaction, proportional 
to the "exchange coupling" J, is approximated here by 
the sum of the squares of the densities of the ± clusters. 

To see the analogy between Eq. (|4.1I) and an Ising sys- 
tem, define 



N, 



iV_ 



N* 



rn 



N* 



V 



N* + N 



(4.2) 



The variable to is analogous to a magnetization; here, it 
is the bond-oricntational order parameter, rj is a measure 
of how close the system is to jamming; it vanishes in the 
dilute limit, No — ► oo, and goes to unity at high density 
where the voids are squeezed out of the system. In terms 
of these variables, Eq. (|4.1j) becomes 



V(m, rj) 
N* 



--l) v -]j V (l + m 2 ). (4.3) 
i] ) 4 

Note that the term proportional to J contains a factor 
rj, implying that ordering becomes weaker with increased 
numbers of voids. Note also that V(m, 77) is proportional 
to N*. Because the partial entropy S also must be pro- 
portional to N* , the latter quantity cancels out of the 
formula for the pressure, and there is no need to include 
it among the relevant internal variables. As a result, we 
can assume that {A} consists of just the two variables m 
and 77. 

Equation (|4.3p provides a formula for the packing frac- 
tion 



(v) 

Ntot y , 



(4.4) 



where Ntot is the fixed total number of particles in the 
system, and (v) is the average volume of a single particle. 
If we measure all volumes, including J, in units such that 



Not (v) = N* 



(4.5) 



1 ) vq - 7 Jr](l + m 2 ). (4.6) 



then we can write 
1 
V 



where v = v* + V\/N*, and V\ is the residual volume 
defined in Eq. (|2.6[) . This scaling implies that both v and 
t>o are dimensionless numbers of the order of unity. The 
total number of particles contained in V\ must be pro- 
portional to Ntot, i-e. V\ ~ Ntot (v)', and therefore, with 
the volume units defined in Eq. (|4.5|) . v ~ Vi/N* ~ 1. 
Similarly, the total volume associated with voids must 
scale like Nq tjq ~ Not i v ) ■ Since we have required Nq to 
scale with TV* via Eq. (|4. 2|) , we again may use Eq. (|4.5I) to 
find that Vq ~ 1. 



C. Entropy 

Much of the physics of this model is contained in the 
choice of the entropy <S, defined in Eq. l|2.5l) . Like V, 
<S must be proportional to the number of statistically 
relevant degrees of freedom, N* . Assume that S can be 
written in the form 



S(m,rj) = <Si(m) +5 2 (ry). 



(4.7) 



The two-state model implies that Si is an Ising-like func- 
tion: 

= ln(2)-i(l + rn)ln(l + m) 



— — (1 — m) In (1 — to). 



(4.8) 



The choice of £2 is more interesting and problematic. 
If we make a lattice-gas approximation in which the Nq 
voids are distributed randomly over N* + Nq sites, we 
find 



S2W 

N* 



In fa)- --1 In(l-ry), 



(4.9) 



which has an ideal gas limit as 77 — > 0, but vanishes very 
weakly as rj —¥ 1. On the other hand, van der Waals 
behavior, withp ~ dS/dr] ~ ( 1 — 77) — 1 , would require that 
£2 ~ ln(l — 77), which is unphysical because it diverges 
as rj — > 1. 

My proposed alternative is 



N* 



- In (77) + 



.4 



1 



(I-'?) 1 



(4.10) 



where the parameter e must be in the range < e < 1. 
Along with the adjustable parameter A, e tunes the 
strength of the jamming behavior between weak and 
strong limits. For small e, £2 approximates lattice gas 
behavior; as e — > 1, it resembles van der Waals. In any 
case, the choice of 5 2 for large 77 must be regarded as 
a phenomenological strategy for data fitting. We know 
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that this model lacks the ingredients for describing the 
way in which the system becomes jammed or crystalline 
at high densities. This is not a first-principles theory of 
such behavior; but it is useful to work with an approxi- 
mate theory in which the jamming limit can be described. 



D. Spatial Variations 

A more general formulation of this theory starts with a 
partition function expressed as a functional integral over 
spatially varying values of m and 77: 



Z(X) = 5 



5r) 



(4.11) 



Relations such as Eq. (|2.10j) can be interpreted as mean- 
field results, obtained by making a saddle-point approx- 
imation in Eg. (14.111) . As usual, write T/N* = f(m, 77), 
and add a square-gradient term in the bond-orientational 
order parameter m: 

^r^/K^ + f (Vm) 2 . (4.12) 

Just as in the magnetic Ising model, the square-gradient 
term has its origins in the pairwise interactions propor- 
tional to J in Eq. (|4~T1) . 

The standard procedure [2(| for dealing with criti- 
cal systems of this kind is to start with the unrenor- 
malized form of of f(m,r]), expand it in powers of m 
to obtain a Landau approximation, and then perform a 
renormalization-group analysis. Beyond m 4 , the higher 
powers in the expansion become irrelevant; and it is easy 
to check that the fluctuations in 77 are non-critical. Thus, 
we know that this procedure produces the correct Ising 
scaling exponents for any set of starting parameters such 
that /(m,7y) has a mean-field critical point. Of course, 
this procedure tells us nothing about the possibility that 
the Ising symmetry might cross over to something else at 
large length scales. 



V. EQUATIONS OF STATE 

We already know from [l[ and Q that the correlation 
exponents for large enough polydispersities A are consis- 
tent with the renormalized Ising values v = 2/d. Thus, 
the most interesting comparisons now are for the equa- 
tions of state, p((j>), for different A's, reported in [2|. 

The mean-field approximation using Eq. (|2.10|) means 
that we look for minima of f{m,rf) in the space of vari- 
ables to and 77. Thus, for variations with respect to m: 



P 



1 



In 



1 



8 J r/m y 1 — m 
and, for variations with respect to 77 



P 



1 



V + (J/4) 77 2 (1 + 771 2 ) 



n - 



(1 - v y 



(5.1) 



(5.2) 
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FIG. 1: (Color online) Pressure as a function of packing frac- 
tion for polydispersties A = 9 % (red open squares) and 13 % 
(blue open triangles). The critical points are indicated by 
the corresponding solid symbols. The blue theoretical curve 
shows the J = limit according to Eqs. (|5.3[ ) and (|5.4[) . The 
pressure is expressed in dimensionless units defined in @]. 



The reference volume v defined in Eq. (|4.6p . and the pa- 
rameters vq, A and J, should depend, at least in a first 
approximation, only on A, and not on 77. With known 
values of these parameters, Eqs. (|4.6p . (|5.ip . and (|5.2[) 
can be solved for p/9, m, and 77 as functions of <f>. 

Note that these equations recover a perfect gas law at 
very small values of 77 and m = 0. Eq. (|5.2p implies that 
p w 9r]/vo; and Eq. (|4.6p implies that <fi w 77/tjo- Thus, as 
expected, p ~ 6(f). 

A second limiting behavior of these equations is espe- 
cially interesting. Large polydispersities A are roughly 
equivalent to high temperatures, which, in the magnetic 
analogy, imply small values of the coupling coefficient J. 
Setting J = in Eqs.(gU) and ([5^]) . we find 



and 



1 



PV 



i-1 



—> 77 



Arj 2 



(5.3) 



(5.4) 



{1-rtY 

If the parameters v, vq, A, and e are independent of A in 
this limit, then the functions p((f>) should collapse onto a 
single curve. I show in Fig.[T]that this is nearly what hap- 
pens for A = 9% and 13%, the latter being the largest 
value for which I have data available. This limiting be- 
havior provides a convenient way to fix several of the 
theoretical parameters. 
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The theoretical curve in Fig. Q] is computed from 
Eqs. (|5.3p and (j5.4p as follows. I arbitrarily have set 
vq = 1, in accord with the argument following Eq. (|4.6l) . 
I also have set 9 = 0.025, which is the value given in [Sj, 
where it is measured in units of the strength of the trun- 
cated, repulsive, Lennard- Jones potential used in those 
simulations, and therefore sets the scale for the pressure 
p. The remaining parameters have been chosen to pro- 
vide an accurate fit to the data - more accurate, I sus- 
pect, than is warranted by either the quality of the theory 
or the precision of the simulations. In this spirit, I find 
e = 0.60, about half way between the lattice gas and van 
der Waals limits, as discussed following Eq. (|4.10l) . Like 
vo and 9, I assume that e is a constant, independent of A. 
According to Eq. (|5.3[) . the maximum packing fraction is 
4>max = 1/v. For this large- A limit, and this choice of e, 



I find 



0.832, and A = 5.6. 



The Ising nature of Eqs. (|4.6[) . (|5.1[) and (|5.2I) appears 
in their reflection symmetry under m — > — m, which is 
spontaneously broken at critical points. We can see this 
behavior at the mean-field level by expanding the log- 
arithm in Eq. (|5.ip to third order in m and solving the 
equation to find 



0. 



where the critical value of rj is 



for J] < r/ c 
for r\ > rj c , 



Vc 



29 

PcJ' 



(5.5) 



(5.6) 



and p c is the critical pressure. This "magnetization" for- 
mula changes substantially under renormalization; the 
mean-field exponent implied by the square root, /3 = 1/2, 
changes to /3 = 1/8 for d = 2 and /3 ^ 0.325 for d = 3. 
This abrupt increase in m for 77 > rj c controls the be- 
havior of the pressure at the onset of ordering via the 
quantity m 2 in the denominator of the right-hand side 
of Eq. (|5.2p . Any meaningful comparison of Eq. (|5.2p with 
the data requires that we use an expression for m 2 that 
is consistent with the renormalized theory. 

To approximate the renormalized behavior, I propose 
to replace the magnetization formula, Eq. (|5.5p . by one 
with the correct Ising exponent: 



M( V ) 



0, 



i)] 



for 77 < rjc 
for 77 > rj c 



(5.7) 



and to use this approximation in Eqs. (|4.6p and ([57 
instead of the mean-field value of m determined by 
Eq. (l5.ip . In the absence of a simple, accurate interpo- 
lation from small to large values of ?7/?7 c — 1, I have used 
Eq. (|5.7p for all 77, and have let pL be an adjustable pa- 
rameter. The only vestige of Eq. (|5.ip in the theory is 
the formula for ?7 C in Eq. (|5.6p . which is obtained from 
Eq. (|5.ip by taking the limit m — > +0. We now may in- 
terpret p c , <fr c , and 77 c to be the renormalized values of 
those quantities; thus they are numbers that we can de- 
duce directly from the data. 
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FIG. 2: (Color online) Pressure as a function of packing frac- 
tion, from bottom to top, for polydispersties A = % (red 
circles), 5% (green triangles), 7% (black circles), and 13% 
(blue triangles). 



The results of these calculations, along with the simu- 
lation data from Q, are shown in Fig.[2]for selected poly- 
dispersities: A = %, 5 %, 7%, and 13 %. For A < 8 %, 
the values of p c and <f> c have been estimated directly from 
the data, because the cusp-like changes in slope at the 
critical points are visible. For larger A's, 4> c has been 
evaluated in [2j by fitting various relaxation-time mea- 
surements to a VFT formula; and the value of (f) c = 0.787 
for A = 9 % also is consistent with the correlation-length 
measurements reported in [![. The critical points for 
A = 9 % and 13 % are indicated respectively by the solid 
square and the solid triangle in Fig. [TJ 

If we know p c and <p c for a given A, then Eqs. (l4.6p 
and (|5.2I) . evaluated at m = and i] = r] c , plus the 
definition of ?7 C in Eq. (|5.6p . provide three constraints on 
the five remaining unknown parameters: ?7 C , v, J, A, and 
fi. A fourth constraint is obtained by estimating 4> m ax 
by setting 77 = 1 in Eq. (|4.6p . but keeping nonzero J; that 
is 



1 



5 [1 + M 2 (1)]. 



(5.8) 



This is an approximate relation, because M(rf) given by 
Eq. (|5.7p is not accurate at ij = 1; but Eq. (|5.8p is a useful 
consistency check on the earlier large-A estimate, espe- 
cially in view of our lack of information about the jammed 
or ordered states in this limit. The only remaining free 
parameter in this analysis is fi. In computing the theoret- 
ical curves shown in Fig 121 I have used \x as the primary 
fitting parameter, and have kept <fimax = 0.832 for all A. 
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The comparisons between the theoretical curves and 
the data in Fig 121 and similar comparisons for other 
values of A not shown here, reveal physically plausi- 
ble trends in the underlying parameters, especially the 
coupling coefficient J. The main trend is that the criti- 
cal points shift to higher values of <j> with increasing A, 
because increasing polydispersity suppresses the order- 
ing transition. The theoretical mechanism that produces 
this effect is the decrease in J shown in Fig. [3] In accord 
with the observations of Tanaka and coworkers, J drops 
abruptly at about A = 8 %; but it does not drop to zero. 
The values of J for A = 9% and 13% have been com- 
puted as above, but using values for 4> c obtained from the 
correlation length and relaxation time analyses instead of 
directly from the p(<j>) curves. This procedure produces 
values of v and A that are consistent with the J = 
analysis, and does not visibly affect agreement with the 
data for p(<p); thus it serves as a self-consistency check. 
The small but nonzero values of J at large A produce 
Ising-like glass transitions. 
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FIG. 3: (Color online) Dimensionless coupling coefficient J as 
a function of percentage polydispersity A 



There are also more subtle effects. The reference vol- 
ume v should be a slowly decreasing function of A be- 
cause, at large polydispersities, smaller particles can fit 
into the spaces between larger particles in ways that do 
not happen in monodisperse systems. The computed v 
does decrease slightly from 1.33 at A = 0% to 1.21 at 
A = 13 %. The parameter A is a measure of the entropy 
and, accordingly, increases with "temperature" A in this 
region, from 4.0 to 5.6, the latter value being the same as 
the one deduced from the J = analysis. The magnetiza- 
tion parameter p, decreases from 1.6 at A = % to about 
1.0 at A = 8%, but becomes impossible to measure at 
the larger A's where the critical transition is not visible 
in the p(4>) curves. In plotting the curves for A = 9 % 
and 13 %, I have simply set /i = 0. 



VI. RELAXATION RATES 

In Kawasaki and Tanaka devote a sizable part of 
their effort to various measurements of the structural re- 
laxation time r a , which they fit to a VFT function to 
determine critical packing fractions <p c . They assume, as 
do I, that this <fr c is the same as the one that governs 
the diverging correlation lengths. The obvious question, 
then, is whether there is a sense in which the VFT for- 
mula is a manifestation of Ising-like symmetry. Here is 
one suggestion about how that may happen. 

Assume that structural rearrangements in glassy ma- 
terials near thermodynamic equilibrium are governed by 
the appearance, disappearance, and internal transitions 
of shear transformation zones (STZ's) [ill, EE HI • The 
STZ's are localized, two-state, flow defects, which must 
contain excess volume (rather than energy) in these hard- 
core systems, and which can undergo shear transitions 
between their two states in the presence of thermal or 
mechanically generated noise. If the characteristic excess 
volume of an STZ is vz, then the equilibrium population 
of such defects is proportional to a Boltzmann factor 
exp(— vz/X), where X = 0/p. To estimate a sponta- 
neous formation rate, multiply this Boltzmann factor by 
an attempt frequency, p(X)/tq, where To is a microscopic 
time determined by the kinetic energies of the particles 
or the thermal fluctuations of the fluid in which they 
are suspended. The dimensionless factor p(X) describes 
glassy slowing down as X decreases. 

Kawasaki and Tanaka 2] show by direct imaging that 
relaxation events occur primarily in disordered regions, 
consistent with the observation of Widmer- Cooper and 
Harrowell 30] that particles undergo rearrangements in 
regions of high "propensity." In the present picture, this 
observation means simply that the STZ formation vol- 
ume vz is smaller in the disordered regions than in the 
ordered ones, so that the STZ's appear most frequently 
in the former. However, the attempt frequency p(X) 
should be a property of the system as a whole, rather 
than being determined just by the local environments of 
a few particles. (See [28J for a discussion of the role of 
STZ's as dynamic heterogeneities and their relation to 
stretched-exponential relaxation. ) 

A scaling argument made in a different context by 
Kirkpatrick et al. @ is useful for estimating a system 
dependent p{X). Consider a correlated volume of size 
Vcorr ~ £ d - Near a critical glass transition, such a vol- 
ume consists of a large number of particles individually 
making small, random displacements. Therefore, it un- 
dergoes Gaussian fluctuations of magnitude proportional 
to the square root of its size; that is, 5V CO rr ~ S, d ^ 2 - The 
frequency of such fluctuations is proportional to 

p(X) ~ e" s v - /x ~ e- « d/2 /^ ~ e- V*" , (6.1) 

where hyperscaling tells us that w = dv /2 = 1 — a/2 = 1 
for both d = 2 and 3. (Here, a denotes the specific heat 
exponent, which is known to be small.) Thus, we recover 
the VFT formula. 
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We can push this argument a bit further by noting that 
it implies 

l0 4-)*7T^-> D= Pc {a^/e, (6.2) 

where a is the linear size of a particle. The bare corre- 
lation length £o) defined in Eq. (|4.I2|) . may be approxi- 
mately the linear size of a binary cluster, and therefore 
ought to be a small multiple of a. We know that p c in- 
creases with A. Thus, the "fragility" parameter D is pre- 
dicted to increase with A - the glass becomes "stronger" 
- in at least qualitative agreement with the increasing 
values of D shown in the inset to Fig. 7 of @. 

I find the argument leading to Eq. (|6.1[) to be too sim- 
ple, especially in view of the fact that the VFT formula 
has not proved to be an accurate representation for wide 
ranges of ex per imental data. In recent applications of 
STZ theory fl3l [29j], my coworkers and I have found it 
better to determine p(X) from simulation data than to 
try to predict it theoretically. Nevertheless, this argu- 
ment may provide yet another element of support for the 
Ising model of glassy behavior. Perhaps the limitations of 
this model, to be summarized briefly in Sec. lVIIl also will 
help us understand the limitations of the VFT formula. 

VII. CONCLUDING REMARKS 

My main hypothesis is that the orientational bias of 
a population of twofold-degenerate "binary clusters" is 
the statistically relevant internal state variable that de- 
scribes a broad class of disordered systems in which the 
constituent particles have a tendency to develop topo- 
logical order. If and when this hypothesis is correct, the 
system exhibits Ising-like behavior; in particular, corre- 
lation lengths associated with the favored topology ap- 
pear to diverge with Ising-like critical exponents. These 
Ising correlations have been observed in several numeri- 
cal simulations, primarlily by Tanaka and coworkers. As 
suggested in Sec. lVIl the observed Vogel-Fulcher-Tamann 
behavior of structural relaxation times may also be evi- 
dence in favor of the two-state hypothesis. 

Tanaka's sequence of two-dimensional, hard-disk tran- 
sitions, visible in the functions p((f>) shown here in Fig. [2] 
indicates a crossover from hexatic to glassy transitions 
with increasing polydispersity A. As discussed at the end 
of Sec. IIII1 however, we know that the theory is wrong 
in the close vicinity of these critical points, because the 
Ising symmetry must change to full circular symmetry 
when the correlations become sufficiently lon g ra nged. 
We know from numerical studies by Jaster [3l| that 
the monodisperse hard-disk system undergoes a KTHNY 
transition with a correlation length that grows exponen- 
tially near <f> c , in contrast to the power-law growth char- 



acteristic of Ising systems. We also know from Anderson 
et al. [32| that, when examined numerically with high 
precision extremely close to <p c , this transition is revealed 
to be weakly first order. There is also an experimental 
study by Han et al. [33(, who see KTHNY behavior for 
two-dimensional colloids at 3 % polydispersity. 

When looked at somewhat less closely, however, the 
two-state Ising theory appears to be remarkably success- 
ful. It correctly predicts a sequence of critical order- 
ing transitions with diverging correlation lengths, even 
for A = 0. The arguments in Sec. IIII1 if correct, make 
the Ising symmetry seem robust; there is no place in the 
two-state picture for a symmetry-breaking analog of a 
magnetic field. The theory also makes roughly credible 
predictions for the ordered states at cj> > (f> c , where the 
correlations again become short ranged, and mean-field 
approximations may regain validity. The small, negative 
values of dp/d(f) in the transition regions might be phys- 
ically realistic indications of the weak, phase-separation 
instability reported in (32j . 

The conjectured validity of mean-field approximations 
in the ordered regime might make it possible for some 
missing ingredients of the theory to be restored within 
the Ising-like formulation. For example, the theory in 
its present state contains no hint of translational order. 
It does not tell us how or where to look for competi- 
tion between glass formation and crystal growth. It re- 
sorts to a phenomenological expression for the entropy, 
in Eq. (|4. 10[) , for computing the pressure at high packing 
fractions, where translational order should be present, at 
least for small A. It says nothing specific about the orien- 
tations of the local topologies, or the possibility of "grain 
boundaries" between topologically oriented regions. It 
contains no information about how particles of different 
sizes are distributed spatially. For example, in a bidis- 
perse system of hard disks, Donev et al. [34J found phase 
separation between the large and small particles in quili- 
brated structures at the highest densities. If the two- 
state, Ising-like model does provide a reasonable starting 
approximation, then it might accomodate some of these 
other physical properties of glassy materials. 
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